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We study spin and charge currents pumped by precessing magnetization of a single ferromagnetic 
layer within _F|/|7V or F\I\F (F-ferromagnet; /-insulator; A^'-normal-metal) multilayers of nanoscale 
thickness attached to two normal metal electrodes with no applied bias voltage between them. Both 
simple one-dimensional model, consisting of a single precessing spin and a potential barrier as the 
"sample," and realistic three-dimensional devices are investigated. In the rotating reference frame, 
where the magnetization appears to be static, these junctions are mapped onto a four-terminal dc 
circuit whose effectively half-metallic ferromagnetic electrodes are biased by the frequency hco/e 
of microwave radiation driving magnetization precession at the ferromagnetic resonance (FMR) 
conditions. We show that pumped spin current in F\I\F junctions, diminished behind the tunnel 
barrier and increased in the opposite direction, is filtered into charge current by the second F layer 
to generate dc pumping voltage of the order of ~ 1 /iV (at FMR frequency ~ 10 GHz) in an open 
circuit. In _F[/|A'^ devices, several orders of magnitude smaller charge current and the corresponding 
dc voltage appear concomitantly with the pumped spin current due to barrier induced asymmetry in 
the transmission coefficients connecting the four electrodes in the rotating frame picture of pumping. 
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I. INTRODUCTION 

The pursuit of "second generation" spintronic devices^ 
has largely been focused on harnessing coherent spin 
states and their dynamics in metals and semiconduc- 
tors. This requires to maintain and control spin orien- 
tations transverse to externally applied or internal mag- 
netic fields. The salient example of phenomena involving 
both coherent spins and their time evolution is the spin- 
transfer torque where spin current of large enough den- 
sity injected into a ferromagnetic layer either switches 
its magnetization from one static configuration to an- 
other or generates a dynamical situation with steady- 
state precessing magnetization.^ In the reciprocal effect, 
termed spin pumping because it occurs in setups without 
applied bias voltage,'^ microwave driven precessing mag- 
netization of a single ferromagnetic layer under the FMR 
conditions emits pure spin current (not accompanied by 
any net charge flux^) into adjacent normal-metal layers. 
In the conventional picture of spin pumping,^ F\N inter- 
face pumps spin current in both directions,^ so that its 
magnitude is determined by the interfacial parameters 
which govern transport of spins that are noncoUinear to 
the magnetization direction at the interface. ^^'^ 

The spin current emitted from the F layer with moving 
magnetization has been observed^* in early experiments 
only indirectly as an enhancement of Gilbert damping® 
of magnetization dynamics in inhomogeneous structures 
due to the presence of i^|A^ interfaces and fast relaxation 
of pumped spins in good "spin sink""^ TV layers which 
ultimately leads to a loss of the angular momentum.^ 
Very recently it has been converted^^ into the conven- 
tionally measurable voltage signals through the inverse 
spin Hall effect (where longitudinal spin current injected 
into a metal with spin-orbit couplings generates trans- 



verse voltage between lateral edges of the sample*). An- 
other electrical scheme is based on Ni\F\N2 multilayers^^ 
where different voltages develop at different F\Ni inter- 
faces due to backflow spin current (driven by the spin 
accumulation in Ni layers built up by directly pumped 
spin current), which is detected by the precessing F layer 
itself. These experiments suggest that spin pumping de- 
vices could be exploited as generators'^'^ of elusive pure 
spin currents,* where spin current injected from F into 
adjacent N layers carries fast precessing spins in giga- 
hertz range of frequencies offering new functionality for 
metal spintronics.® They can also be used to probe im- 
portant aspects of spin dynamics in thin F layers. 

Unlike spin-transfer torque that has been demon- 
strated in F\I\F magnetic tunnel junctions (MTJ),^ 
it has been considered that low transparent interfaces 
would completely screen the interfacial spin pumping ef- 
fect (as observed in some experiments^*), unless the tun- 
nel barrier has non-trivial magnetic properties. Thus, 
recent surprising measurements"'^^' of large voltage sig- 
nals of the order of ^ 1 fiV (at FMR frequencies / — 2 
GHz and precession cone angles 9 ~ 10°) in microwave 
driven F\I\N and F\I\F tunnel junctions, as opposed to 
~ 10 nV pumping signals"'"'^ in Ni\F\N2 multilayers, have 
attracted considerable theoretical attention. ^®"^° Never- 
theless, the puzzle of unexpectedly large magnitude of the 
observed dc pumping voltages persists: (i) the scatter- 
ing^® approach to transport of noninteracting quasipar- 
ticles through defect-free epitaxial F\I\F MTJ finds ^ 
1 nV signals at FMR frequency f — 2 GHz and preces- 
sion cone angle 9 = 10°; (ii) the tunneling Hamiltonian 
approach"'^ for clean F\I\F MTJ and the same / and 
9 parameters sets the maximum dc pumping voltage at 
^ 0.01 /iV in parallel and ^ 1 /iV in antiparallel con- 
figuration of two F electrodes; and {Hi) the tunneling 
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FIG. 1: (color online), (a) The ID model of spin pumping 
where the sample consisting of two sites, one hosting the single 
spin rotating with frequency ui and the other one hosting the 
potential barrier of height ei, is attached to two semi-infinite 
tight-binding chains (7 is the hopping parameter) playing the 
role of electrodes with no applied bias voltage between them. 
In the rotating reference frame the spin is static and the de- 
vice (a) is mapped into the four-terminal dc circuit in panel 
(b) whose electrodes have electrochemical potential shifted 
by ±hu}/2 with respect to the equilibrium Fermi level Ef of 
unbiased electrodes in the laboratory frame. 



Hamiltonian approach combined with semiclassical mod- 
eling of the interplay of spin diffusion and self-consistent 
screening around interfaces in F\I\F and _F|J|A'' junc- 
tions involves too many unknown phcnomcnological pa- 
rameters, thereby offering only a wide range of possible 
pumping voltages for both of these junctions.^" 

Here we address the problem of spin pumping and in- 
duced voltages by high frequency magnetization dynam- 
ics in F\I\F and i^|/|Af junctions within the framework of 
nonequilibrium Green function (NEGF) formalism. ^^'^^ 
We note that NEGFs have been utilized before to study 
spin^^'^"* and charge^^ pumping by time-dependent fields 
acting on finite-size paramagnetic devices attached to 
electrodes held at the same electrochemical potential. 
Since NEGF formalism takes as an input a microscopic 
Hamiltonian, it makes it possible to include, in a con- 
trolled fashion, the full geometry^^ of experimental de- 
vices (such as the finite thickness of F, /, and N lay- 
ers, down to few atomic monolayers, which play an im- 
portant role in the magnetoresistance^^ and spin-transfer 
torque'^^ of crystalline MT Js) , the properties of the insu- 
lating barrier (including disorder effects), as well as the 
interactions responsible for spin-flip processes in F. The 
NEGF formalism also makes it easy to take into account 
ah initio input^^"'^" on the F\I interface electronic and 
magnetic structure and the self-consistently developed 
nonequilibrium spin and charge distributions around it. 

Furthermore, NEGF approach yields a remarkably 
transparent physical picture of pumping in ferromagnetic 
multilayered systems. For example, in the simplest model 



of pumping, generated by a single spin precessing with 
frequency oj in Fig. 1(a), the NEGF rotated into the rest 
frame of the spin maps the original laboratory-frame de- 
vice onto a dc circuit in Fig. 1(b). The central sample of 
this circuit, which contains time-independent spin inter- 
actions, is attached to four electrodes that allow only 
one spin species to propagate through them and are, 
therefore, labeled by L-left, i?-right, spin-|, and spin- 
l- These four electrodes are biased by the voltage huj/e, 
so that spin- J, electrons flow from electrodes at higher 
electrochemical potential and precess inside the sample 
due to spin-dependent interactions to be able to enter 
into electrodes at a lower electrochemical potential as 
spin-t states. Thus, this picture reduces the quantita- 
tive analysis of spin and charge pumping by precessing 
spins to multiterminal Landauer-Biittiker-type formulas 
for spin-resolved charge currents as encountered in, e.g., 
the mesoscopic spin Hall effect. ^^'^^ 

The paper is organized as follows. In Sec. II we exploit 
the physical picture of pumping provided by Fig. 1(b) 
to analyze local spin and charge currents flowing away 
from the single precessing spin toward the neighboring 
sites along the tight- binding chain in one-dimension (ID). 
This framework is extended to total pumped currents and 
associated voltages in three-dimensional (3D) multilay- 
ered structures, such as F\N\F, F\I\F, and F\I\N, in 
Sec. HI. We conclude in Sec. IV. Our principal results — 
pumped spin and charge currents in ID model and volt- 
age signals in F\I\F and F|/|iV junctions — are shown in 
Figs. 2 and 5, respectively. 

II. NEGF APPROACH TO SPIN AND CHARGE 
PUMPING BY A SINGLE PRECESSING SPIN IN 
ONE DIMENSION 

The toy ID model in Fig. 1(a) encodes most of the 
essential physics of pumping by precessing spins while 
making it possible to obtain analytical solution for the 
magnitude of pumped currents. For simplicity, we start 
from the often employed in spin-transfer torque'^^ and 
spin pumping^^ studies Stoner-type Hamiltonian'^'^ 

r,cr,cr' ^ 

-7 ^LcrV, (1) 

(rr')cr 

in the local orbital basis suited for NEGF calcula- 
tions. ■^^'^^ Its time dependence stems from the unit vec- 
tor m.(t) along the local magnetization direction, which 
is assumed to be spatially uniform and steadily precess- 
ing around the z axis with a constant cone angle 6. The 
operators cj^. (cra) create (annihilate) electron with spin 
cr at site r, and 7 is the nearest neighbor hopping. The 
coupling of itinerant electrons to collective magnetic dy- 
namics is described through the material-dependent ex- 
change potential Ar, where & = [a^, o'y,az) is the vector 
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FIG. 2: The (a) Sz-spin and (b) charge currents pumped by a 
single precessing spin as a function of the potential barrier on 
the second site of the sample in ID model shown in Fig. 1(a). 
The parameters of the model are: / = cj/27r = 20 GHz; 
= 10°, A/Ef = 0.85, and electrons in the macroscopic 
reservoirs to which the electrodes are attached have the Fermi 
energy Ef = 2'y. 



of the Pauli matrices and af^ denotes the Pauli matrix 
elements. The on-site potential Sr accounts for the pres- 
ence of the barrier [such as = £i on the second site of 
the sample in Fig. 1(a)], disorder, external electric field, 
and it can also be used to shift the band bottom conve- 
niently. The sample is attached to two semi-infinite ideal 
(spin and charge interaction free) electrodes, which ter- 
minate in macroscopic reservoirs held at the same elec- 
trochemical potential /ip = Ep where Ep is the Fermi 
energy. 

The fundamental objects^^ of the NEGF formalism are 
the retarded 

G:':/{t,t') = -ie(t-o({cr.(t),4v'(i')}>, (2) 

and the lesser 

G<r\t^t') = \{^r'At')Cr.{t)), (3) 

Green functions ((...) denotes the nonequilibrium sta- 
tistical average^-'^) which describe the density of avail- 
able quantum states and how electrons occupy those 



states, respectively. Since nonequilibrium problems are 
not time-translation invariant, these Green functions de- 
pend on two time variables separately. However, the 
cumbersome double time dependence of NEGF in gen- 
eral pumping problems^^ can be eliminated^^ for the 
special case of time-dependent potential caused by pre- 
cessing magnetization using the compensating rotation^'' 
of the system described by the unitary transformation 



U 



(for magnetization precessing countcrclock- 
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wise). Thus, the Hamiltonian in the rotating frame 

^,ot = UHuhm^ - if^U-U^ = Aab(O) - —a,, (4) 

is time-independent. The term hiuaz/2, which appears 
uniformly in the Hamiltonian of the sample or N elec- 
trodes, will spin-split the bands of the N electrodes. This 
yields a rotating-frame picture of pumping based on the 
four-terminal device in Fig. 1(b). 

The device in Fig. 1(b) guides us in setting up the 
NEGF equations for the description of currents flowing 
between its four electrodes, labeled by p, cr (p ~ L, R 
and a =t, i), which are biased by the voltage hu)/e. The 
electrodes behave effectively as the half-metallic ferro- 
magnets, emitting or absorbing only one spin species. 
The rotating frame Green functions 



and 



G<{E) = G'^iE)^:<{E)G''{E), 



(5) 



(6) 



depend on t = t — t', or energy E after the time differ- 
ence T is Fourier transformed. Here the advanced Green 
function is G'^{E) = [G^{E)]'', and Hrot is the matrix 
representing -ffiot in the local-orbital basis. The retarded 
self-energy matrix S'' (E) = J2p a ^p'^ i^) ''U™ of 

self-energies introduced by the interaction with the leads 
which determine escape rates of spin-cr electron into the 
electrodes p, a in Fig. 1(b). 

For interacting systems Y1'^{E) would also con- 
tain electron-electron and electron-phonon contributions, 
while for noninteracting systems, described by Hamilto- 
nian (4), the lesser self-energy is expressed in terms of 

i:;'-{E) as 



^^<{E)=J2^fpiE)K{E). 



The level broadening matrix 



T;{E) = -2Im'El[E + a 



(7) 



(8) 



is obtained from the usual self-energy matrices^^ ^pi^) 
of semi-infinite leads in the laboratory frame with their 

energy argument being shifted by ahw/2 to take into ac- 
count the "bias voltage" in accord with Fig. 1(b). The 
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distribution function of electrons in the four electrodes 
of the rotating frame dc circuit is given by 



P AP 



fpiE) 



1 



exp[{E -Ef + ahuj/2)/kT] + l' 



(9) 



where a = + for spin-| and a = — for spin-j. Since the 
device is not biased in the laboratory frame, the shifted 
Fermi function in Eq. (9) is uniquely specified by the 
polarization t or | of the electrode, so that we remove 
the lead label p from it in the equations below. 

The basic transport quantity for the rotating-frame dc 
circuit is the spin-rcsolved bond charge current^^ carry- 
ing spin-(T electrons from site r to neighboring site r' 



This gives spin 



and charge 



1 , — f +7^-^ 



(10) 



(11) 



(12) 



bond currents flowing between neighboring sites. Equa- 
tion (10) can be evaluated analytically for ID model 
assuming that for small enough huj Ep we can use 
f^{E) - f^ {E) = hujS{E — Ep) at zero temperature. 
Such "adiabatic approximation"^^ is analogous to linear- 
response limit in conventional transport calculations for 
devices biased by small voltage difference. 

The ^z-component of the bond spin current between 
the precessing site and its nearest neighbor in the sample, 
as illustrated in Fig. 1(a), is given by 

js. , . 2.7'A2(ImEiD)^ 



X [4 (72 + ej) + 4|S]iDp - 8e/Re Em] , (13) 

where the terms O [{hujY] are neglected. Here Sid = 
{Ep - 27 - y/{Ep - 27)2 - 472)72 is the self-energy of 
ID semi-infinite lead (i.e., tight-binding chain) and R ~ 
(4E2^ - A2)(I]iD - e,)2/4 + (4£,EiD - AJ:l^hy2 + 7^. 
In deriving Eq. (13), we assume uniform band bottom, so 
that £1 ^ Si - (Acos6')/2 and Ep ^ Ep + (Acos6')/2 
is used to plot Fig. 1. 

The expression in Eq. (13) reproduces all major fea- 
tures of the scattering approach'^ to adiabatic {hjuo ^ A) 
regime of spin pumping by interface in 3D multilay- 
ers: (j) the pure spin current carrying Sz spins is propor- 
tional to hui and sin^ 6\ {it) Sz component J^^, of the spin 
current tensor is time independent in both rotating and 
laboratory frames; and (iii) Sx- and 5y-components of 
the pumped spin current oscillate harmonically with time 
in the laboratory frame. Moreover, when potential bar- 
rier ej is introduced into the sample, we find in Fig. 2(a) 




FIG. 3: (color online). The magnetic tunnel junction with 
precessing magnetization in the left F layer is modeled on 
a simple-cubic tight-binding lattice. The thicknesses of the 
ferromagnetic {F, F') and thin insulating (J) layers are mea- 
sured using the number of atomic monolayers dp, dpi, and 
di, respectively. The P and AP configurations of MTJ corre- 
spond to the magnetization of the right F layer being parallel 
or antiparallel to the z-axis around which spatially uniform 
magnetization of the left F layer steadily precesses with a 
constant cone angle 6. 



that spin current on the right decays with increasing £/ 
while pumped spin current flowing on the left increases 
to about twice the value of the sum 
the left and right spin currents pumped symmetrically in 
the absence of the barrier. This effect can clarify the ori- 
gin of possible Gilbert damping enhancement in realistic 
MTJ devices consisting of iV|F|/|i^| multilayers, rather 
than infinite F electrodes, where angular-momentum loss 
develops due to increasing spin pumping into the left N 
electrode even when the insulating barrier / suppresses 
spin pumping on the right side of the junction. 

Figure 2(b) demonstrates that non-zero potential ej ^ 
also leads to concomitant pumping of a tiny charge 
current into the right electrode, which is several order of 
magnitude smaller than the pumped spin current. Un- 
like pumping of spins which is linear in frequency, such 



pumped charge current scales as 



We discuss its 



origin in Sec. Ill by analyzing total charge current in the 
N terminals expressed in terms of the transmission coef- 
ficients between the fully spin-polarized electrodes in the 
rotating frame. 



III. NEGF APPROACH TO SPIN AND 
CHARGE PUMPING IN 3D F\N\F, F\I\F, AND 
F\I\N MULTILAYERS 

We extend this analysis to a 3D MTJ shown in Fig. 3 
which consists^^^ Qf infinite planes of F, N, and / materi- 
als modeled on a simple-cubic tight-binding lattice with 
single s-orbital per site using Hamiltonian (1). The effec- 
tive dc circuit [Fig. 1(b)] in the rotating frame makes it 
easy to write the expression for the total charge current 
in the left and right N electrodes. For example, spin-| 
electrons can flow from ^ lead at higher electrochemical 
potential into ]^ lead at the lower electrochemical poten- 
tial. They enter ]j lead as spin-| electrons with probabil- 
ity determined by the precession inside the sample since 
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-ffrot contains terms proportional to a.j. for which the in- 
jected spin states ||) from lead (polarized along the z 
axis) are not the eigenstates. 

Since pumped charge current is necessarily conserved, 
as exemplified by Fig. 2(b), we arbitrarily select the right 
electrode (current flowing into the electrode is assumed 
to be positive) to find its explicit expression in terms of 
the multiterminal Landauer-Biittiker formulas^^ for spin- 
resolved quantum transport: 

oo 

/ = ^ JdE{Tj,i[f^iE)-fiE)] 

— OO 

-Tl\,[fl{E)-f{E)]}. (14) 

Here the transmission coefficients T^^, determine the 
probability for a' electrons injected through lead p' to 
emerge in electrode p as spin-a electrons. They can be 
computed from the NEGF-based formula^^ 

which is written here in the spin-resolved form. The block 
Gp^i"^ of the retarded Green hmction matrix consists 
of those matrix elements which connect the layer of the 
sample attached to lead p' to the layer of the sample 
attached to lead p. 

In general, the spin current is not conserved, as illus- 
trated by Fig. 2(a), and we choose to compute it in the 
left A'' electrode: 

oo 

It = I jdE{Tl^^[f^{E)-f\E)] 

— OO 

+Tli[fHE) - fHE)] 

+ 2Tli[f^{E)-f{E)]}. (16) 

The expressions Eq. (14) and Eq. (16) for total currents 
are equivalent to the sum of all bond charge Jrr' or bond 
spin J^, currents, respectively, where summation is per- 
formed over the pairs of sites within the electrode at a 
chosen cross section. By the same token, the analyt- 
ical expression Eq. (13) for the spin current is already 
equivalent to the result obtained from Eq. (16) since no 
summation is necessary for the cross section consisting of 
a single site. 

The pumped charge current in multilayers with the 
second analyzing F layer originates from spin filtering 
by the static magnetization of the analyzing F layer of 
current pumped toward the right. That is, we find I^' — 
0, I^" 7^ and / 7^ in such systems. In junctions with a 
single precessing F layer the; pumped spin current is pure 
if jf * = ^ and / = 0. The possibility of nonzero 
pumped charge current even in junctions with only one 
F layer whose magnetization is precessing, as exemplified 
by Fig. 2(b) and F| J| junctions in general [see Figs. 5(c) 



and 5(d)], is explained by Eq. (14) as the consequence of 
the asymmetry in transmission coefficients 1^^/,— I^jj 7^ 
when arbitrary potential £/ ^ is introduced in one of 
the layers. 

The transmission coeflScients can also explain the un- 
expectedly large enhancement of /f ' or J^^^^ in Fig. 2. 
As the barrier height £/ increases, and tJj|^ dimin- 
ish to very small value while 2T2^ increases to about four 
times its value at £/ = due to quantum interferences ef- 
fects on the left side of the device (quantum interferences 
were also found to enhance pumped spin current when co- 
herent backscattering from disorder occurs in finite-size 
conductors at paramagnetic resonance^^). At e/ ~ 0, 
T^^ = T^;^ = t]^j^ so that spin currents of the same 
magnitude are pumped in both directions symmetrically. 

The pumped charge current is translated into dc volt- 
age in open circuits via 

^ump = 'q(^i (■'■'^) 

where G(^) is the conductance of F\1\F (or i^|/|A'' when 
the second F layers is removed) junction sketched in 
Fig. 3 whose first F layer has its static magnetization 
tilted by an angle away from the z ax;is and the linear 
response bias voltage is applied between the electrodes 
in the laboratory frame. The quantity G{d) = 2e^Tf{i//i 
can also be computed via the standard NEGF formula^ ^ 
as in Eq. (15) but for total Tjil, rather than spin- 
resolved, transmission coefficient expressed in terms of 
the retarded Green function and self-energies in the lab- 
oratory frame. 

The largest voltage signal of spin pumping is ex- 
pected in high quality epitaxial Fe|MgO|Fe tunnel jmic- 
tions.^® To mimic their huge tunneling magnetoresis- 
tance (TMR), while using the simple single-orbital tight- 
binding Hamiltonian (1), we adopt the same parameters 
employed in Ref. 18: Ep = 4.5 eV, A/i;^. = 0.85, 7 = 1.0 
eV, and the barrier height measured relative to the Fermi 
energy Ub = {ei - Ep) is Ub/Ep = 0.25. The band bot- 
tom is aligned across all layers of the junction with the 
bottom of the band for majority spins in F (similarly 
to Ref. 18). The "optimistic" TMR ratio for this junc- 
tion with di = 5 monolayers of the insulating material is 
TMR = {Rap - Rp)/Rp ^ 3900%, which is close to ab 
initio computed zero-bias TMR ~ 3700% for defect-free 
Fe|MgO|Fe MTJ containing five MgO layers.^° 

In the coherent limit of tunneling, '""^ applicable to ideal 
crystalline structures without any defect scattering, the 
in- plane wave vector k|| = {ky,kz) is conserved and all 
NEGF quantities depend on it. This requires to inte- 
grate T^^,'{E,k\\) in Eq. (14) and (16) over the two- 
dimensional (2D) Brillouin zone (BZ). Thus, in the adia- 
batic limit and at zero temperature we use the following 
formulas to obtain the charge current 

I=^^j<^\\{T'A-Tlil (18) 

BZ 
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FIG. 4: (color online). The dc pumping voltage in _F'|A'^|_F 
multilayers attached to two semi-infinite N electrodes as the 
function of the thickness of F layer whose magnetization is 
precessing with cone angle 9 = 10° at frequency f = uj/2n = 
20 GHz. The parameters describing the multilayer are Ep = 
4.5 eV, A/Ef = 0.85 (in both F layers), and 7 = 1.0 eV. 



and the spin current in the left lead 
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FIG. 5: (Color online). The dc pumping voltage in [(a) and 
(b)] F\I\F and [(c) and (d)] F\I\N multilayers attached to 
two semi-infinite A'' electrodes as the function of the barrier 
height Ub (measured relative to the Fermi energy). The mag- 
netization of the left F layer is precessing with cone angle 
6 = 10° at frequency / = tj/27r — 20 GHz. The parameters 
describing the multilayer are Ef ~ 4.5 eV, A/Ef = 0.85 (in 
both F layers for F\I\F junction), and 7 = 1.0 eV. 



pumped by magnetization precessing at frequency oj. 

The computational algorithm for this integration can 
be substantially accelerated by transforming the 2D pla- 
nar momentum integral into a single integral over the 
in-plane kinetic energy 



7T / a 7r/a 



-TT /a—TT I a 



xT;; {E,ey._), (20) 



where we utilize the two-dimensional density of states 
P2D{£yz) for a square lattice and the fact that T^T de- 
pends on k|| through the in-plane kinetic energy Syz- In 
the case of nearest-neighbor hopping on a square lat- 
tice, the kinetic energy within a monolayer is given by 
£yz — 4:'-f — 2j[cos{kya) + cos(fcza)], where a is the lattice 
spacing. The effect of the in-plane kinetic energy is equiv- 
alent to an increase in the on-site potential £r ^ £r+£yz- 
To provide reference values for understanding the mag- 
nitude of pumping voltages in tunnel junctions, as well as 
to connect our theory to a "standard model" of interfacial 
spin pumping provided by the scattering theory, '^'^■^ we 
first compute the dc voltage Vpump generated in F\N\F 
multilayers. The chosen cone angle 9 = 10° and FMR 
frequency / ~ lu/2tt — 20 GHz are within the range of 
typical values encountered in experiments^^ where the 
results in Figs. 4 and 5 can easily be rescaled for other 
values of these two parameters using the general oc sin^ 9 
and oc lo dependence in Eq. (13). Figure 4 demonstrates 
that pumping involves only a thin layer of F material 



around the F\N interface. However, while in the scatter- 
ing theory'^ adiabatic pumping develops over the atom- 
istically short ferromagnetic coherence length ~ hvp/^, 
which in our junction is Hvp/^ — a, we find that pump- 
ing in Fig. 5 involves about five monolayers of the ferro- 
magnetic material. Here we assume that the magnitude 
of pumped current generated on this length scale is not 
affected^^ by spin-relaxation processes [not included in 
Hamiltonian (1)] that typically occur on a much longer 
length scale. '^^ The pumped voltages in both P (parallel) 
and AP (antiparallel) configuration are below the maxi- 
mum"^'^* expected voltage Vpump < fi^ « 83 fiV (for the 
explanation of P and AP junction setups in the context 
of pumping by precessing magnetization, see Fig. 3). 

The dc pumping voltage for tunnel junctions is shown 
in Fig. 5. Although the presence of the potential barrier 
within / layer of junction increases the resistance 

of the junction in Eq. (17), the pumped charge current 
decreases faster so that Vpump decreases with increasing 
barrier height Ub- In contrast to the scattering result of 
Ref. 18 where Vpump increases with increasing Ub for all 
thicknesses of the / layer, we find in Fig. 5(b) such in- 
crease only if the I layer consists of a single monolayer. 
The large difference between Vp^j^jp and Vp'^^p configura- 
tion stems from huge TMR ratio for this junction while 
the magnitude of pumped charge current remains virtu- 
ally the same for both P and AP configurations. 

These results are quite close to 1 /iV for Vp^^p — 
^p^imp observed in MTJs with AI2O3 barriers, and an 
order of magnitude larger voltages in MTJs with MgO 
barriers. The experimentally observed^^ change in sign 
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of Vpump depending on the type of the barrier (AI2O3 vs. 
MgO) or its thickness might be related to a difference 
in sign between Figs. 5(a) and 5(b). We also compute 
the pumped spin currents in the left ifj' = 0.012e//i and 
right = electrodes, which do not depend on £/ or 
Ub in the range shown in Fig. 5. 

Analogously to pumped charge current of ID model 
in Fig. 2(b), we find nonzero charge current and corre- 
sponding dc pumping voltage in _F|/|A^ junctions shown 
in Figs. 5(c) and 5(d). Nevertheless, Vpump of the order of 
~ 10 pV are way to small to explain recent experiments 
on i^|/|A'' junctions^® where T4>ump — 1 A^V is measured 
at frequencies of the applied rf field in the range / = 2-3 
GHz and the precession cone angle 6 = 10°-17° tuned by 
the microwave input power. 

IV. CONCLUDING REMARKS 

In conclusion, we have demonstrated that pumping of 

spin and charge currents by the processing magnetization 
of a ferromagnetic layer within various multilayer setups 
consisting of F, N, and / layers of nanoscale thicknesses 
can be understood within the framework of NEGF ro- 
tated into the frame moving with the magnetization as 
a simple four-terminal dc circuit problem, as illustrated 
by Fig. 1(b). The four leads of this circuit are labeled as: 

^, ]j, and (i.e., they act as half-metallic ferromag- 
netic electrodes). They arc biased by the voltage differ- 
ence hui/e effectively emerging between the electrodes of 
opposite polarization. Our formalism provides a trans- 
parent physical picture of how: (i) single precessing spin 
pumps pure spin current symmetrically (in the absence 
of any barriers) toward the left and the right in ID; {ii) 
pumped spin currents are suppressed by the tunnel bar- 
rier in one direction and enhanced in the opposite direc- 
tion beyond naive sum of currents before the introduction 
of the barrier; (iii) pumped spin currents develop over 
few monolayers of F material in 3D junctions; and (iv) 
pumped spin currents become filtered by the second F 
layer with static magnetization which converts them into 
charge current and the corresponding dc pumping volt- 
age in open circuits. Our physical picture of spin and 
charge pumping in MTJs with time-dependent magneti- 
zation suggests that these setups can serve as a sensitive 
probe of MTJ parameters, such as the properties of the 
tunnel barrier and damping parameters. 

The pumping voltages in A'^ji^j J|F|A'' tunnel junctions 
of the order of ^ 1 /xV at FMR frequencies 10 GHz 
could explain some of the recent measurements of large 
voltage signals in microwave driven MTJs under the 
FMR conditions. They are much larger than ~ 10 nV 
signal (at FMR frequencies 10 GHz) recently pre- 



dicted by the scattering theory^^ for MTJs with simi- 
lar TMR, but whose infinite F electrodes are assumed 
to have strong spin-flip scattering leading to a vanishing 
spin accumulation in F}^ The spin- flip scattering can 
easily be introduced in Hamiltonian (1) via spin-orbit 
(SO) coupling terms'^^ whose strength is tuned to match 
experimental values for spin-diffusion length.^^ Neverthe- 
less, here we use simpler Hamiltonian following assump- 
tions similar to Ref. 19 typical spin-relaxation lengths'^^ 
are much longer than the length scale (illustrated by 
Fig. 4) over which pumping develops so that it does not 
affect the strength of pumped currents. On the other 
hand, computation of realistic patterns of spin accumu- 
lation^° throughout the device requires to consider bal- 
ance between transport and relaxation processes. Also, 
the NEGF formalisms developed here, with spin-diffusion 
length vs selected layer thickness tuned via microscopic 
SO scattering terms, can tackle complicated spin pump- 
ing multilayer setups involving I layers where conven- 
tional approaches are not applicable (because of spin ac- 
cumulation not being well defined in an insulator). 

Although we do find nonzero charge current in 
A'^|i^|/| A" multilayers when potential barrier is introduced 
in the device through the / layer, the voltage signal ~ 10 
pV is several orders of magnitude smaller than ~ 1 /xV 
observed in experiments on such devices. Also, this 
charge ciurent is proportional to {hu})'^, rather than fiui 
for spin and charge currents in N\F\I\F\N junctions or 
experimentally observed dc voltage signal in N\F\I\N 
junctions. -"^^ Its origin is in asymmetry of transmission 
coefficients connecting the four electrodes of the dc cir- 
cuit in the rotating reference frame. While the magni- 
tude of measured voltages remains a puzzle for a vari- 
ety of approaches'^^"^" utilized very recently to address 
some of the aspects of the experiment in Ref. 16, we be- 
lieve that combining NEGF approach to spin pumping 
outlined here with the density functional theory (DFT) 
to take into account nonequilibrium self-consistent spin 
and charge densities (akin to NEGF-DFT approach^^ to 
spin-transfer torque in spin valves and MTJs) could be 
capable of addressing this problem. 
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